A&A manuscript no. 

(will be inserted by hand later ) 



Your thesaurus codes are: 

08 (19.42.1 HD 100546, SAO 206462, MWC 480; 19.25.1; 04.01.1; 05.01.1; 

19.94.1) 



ASTRONOMY 

AND 
ASTROPHYSICS 

February 1, 2008 



HST/NICMOS2 coronagraphic observations of the 
circumstellar environment of three old PMS stars: 
HD 100546, SAO 206462 and MWC 480. * 



O' 

o: 
o. 

(N ■ 

CW 
<D ■ 
00 

OS 



> 

O 
O 
O 

Oh' 
i 

O 



J.C. Augereau 1 , A.M. Lagrange 1 , D. Mouillet 1 , and F. Menard 1 ' 2 

Laboratoire d'Astrophysique de l'Observatoire de Grenoble, Universite J. Fourier / CNRS, B.P. 53, F-38041 Grenoble Cedex 
9, France 
2 Canada-France-Hawaii Telescope Corporation, PO Box 1597, Kamuela, HI 96743, USA 

Received ; accepted 



* 



Abstract. The close environment of four old Pre-Main 
Sequence stars has been observed thanks to the coron- 
agraphic mode of the HST/NICM0S2 camera at A = 
1.6 fiin. In the course of this program, the detection of a 
circumstella r annulus around HD 141 569 has already been 
presented in Augereau ct al. (1999b ). 

In this paper, we report the detection of an elliptical 
structure around the Herbig Be star HD 100546 extend- 
ing from the very close edge of the coronagraphic mask 



(~50Aj U) to 050 080 AU (3.5 3.8") from the star. Tin 



axis rati 



10 juu jouau (j.u j.M'j ii'om tnc star, 
gives a disk inclination of 51° _L 3 g to the linc- 



1. Introduction 

Dust is present around a large fra ction of Main Sequenc e 
(MS) stars, as shown by I RAS ([Aumann ct al. 1984| ). 
Moreover, spectroscopic (e.g. Queloz et al. 2000 ) and pho- 
tometric ( Charbonncau ct al. 2000 ) detections indicate 
that planetary systems are quite common. In some cases, 
the dust lifetime is so short compared to the star's age, 
that it has to be replenished, probably through collisions 
between planetesimals or through evaporation (Backman 



Paresce 1993). The best example so far among these 



of-sight and a position angle of 161°±5°, measured east of 
north. At 50 AU, the disk has a surface brightness between 



second generation disks is the one surrounding j3 Picto- 



ris aged more than 20Myr (Barrado y Navascues et al 



10.5 a rid 11 mag/arcscc", then follows a — 2. 92 _L 0.04 radial 
power law up to 250-270 AU and finally falls as r - 5 - 5±0 - 2 . 
The inferred optical thickness suggests that the disk is 
at least marginally optically thick inside 80 AU and opti- 
cally thin further out. Combined with anisotropic scatter- 
ing properties, this could explain the shape of a brightness 
asymm etry observed along the minor axi3 of the di3k. This 



1999). Recently, submillimeter images resolved the dust 



emission around a few isolated MS stars with different 



spectral types and older than 0.1 Gyr (Holland ct al. 1998 
preaves et al. 1998) 



Circumstellar disks are also present around young ob- 



jects, as clearly demonstrated by HST (Padgett et al. 1999 



for instance), and by ground-based observations (Dutrey 



asym metry nccd3 to be confirmed. 

The circumstellar disks around SAO 206462 and 



et al. 1996). The latter detected circumstellar disks around 



a sample of CTTs in the Taurus cloud, but not yet around 



MWC 480 are uuL resulved, leading Lu euiisLiahiLs on Lire 



post T Tauri stars. This is in general agreement with Zuck 



dust di^liibuliun. A light binary system sepaiaLed by uiily 
0.32" J.0.04" is nevertheless deleeled in lire eluse vieinily 
of SAO 206462. 



erman & Becklin (1993) and more recently Habing et al 
(1999[) who have shown the rapid decrease of the circum- 



stellar dust mass as the stars evolve to the Main Sequence 



(see also Holland et al. 1998). Gas infall is also observed 
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Bas ed on observations with the NASA/ESA Hubble Space 



towards a few Herbig Ae/Be (hereafter HAeBe) stars (e.g. 
prady et al. 1996 ; j.e Winter et al. 1999| ), tentatively inter- 
preted as the result of comet evaporation. In this scheme, 
this could indicate that planctesimal formation in the neb- 
ula is quite common and rapid, as theoretically expected. 
Very few disks have been detected around HAeBe and 
more evolved Pre Main-Sequence (PMS) stars (HD 163296 
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and MWC 480, pVIannings fc Sargent 1997a| [Mannings et 
al. 1997b| ). A recent study of HR4796, an 8Myr star, has 



evidenced that the disk is at least partly second genera- 
t ion and that meter-size bodies should already be present 
( Augereau et al. 1999a ). Similar conclusions co uld apply 
for the HD 141569 disk QAugcrcau et al. 1999b| ). In both 
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cases, the analysis relies on detailed studies of the resolved 
images and spectral energy distributions (SEDs). From a 
general point view, the nature of transient disks, the role 
of replenishing planetcsimals if any, and the most relevant 
physical processes in these environments are still to be 
clarified. 

The general evolutionary scenario, to be tested and 
precised is then that : 

1) the protoplanetary disks around young, embedded 
stars clear out before the star reaches the main sequence 
through dynamical activity and interaction with stellar 
radiation and that 

2) new optically thin circumstellar disks may be sustained 
around older MS stars by destruction processes among 
planetesimals (collisions, evaporation), formed during the 
Pre-Main Sequence (PMS) phase. 

The transition between those two phases, when the origi- 
nal young disk is eroded, is critical to validate the above 
scenario. Such a study can be performed on disks around 
evolved PMS stars. Critical points for this study are the 
analysis of opacity to stellar radiation and collisional pro- 
cesses time-scales. Both of them require detailed informa- 
tion on the spatial distribution and amount of circumstel- 
lar dust. 

During HST/Cycle 7, we probe in coronagraphic mode 
the environment of four evolved PMS stars with the NIC- 
MOS2 camera. As part of this program, we reported the 
positive detection of an extend ed (> 400 AU in radius ) 
structure around HD 141569 in [Augereau et al. (1999b ). 
The detection was confirmed by Weinberger et al. (199E ) . 
We address in the present paper our results for three other 
sources of large interest among old PMS stars : HD 100546, 
SAO 206462 and MWC 480. After a brief description of the 
observing strategy, journal of observation and reduction 
procedure (section g) , we present in section || the results 
for HD 100546. An extended structure is detected and in- 
terpreted in terms of circumstellar dust arranged within 
an inclined disk. Detection limits are derived from unre- 
solved structures around SAO 206462 and MWC 480 (sec- 
tions and ra) and in each case we discuss the implications 
on the dust distribution. We also detail the close vicinity 
of SAO 206462 which exhibits a tight binary system. 



ized by a large IRAS excess, together with an age estima- 
tion (via photometry or spectroscopy) of about 10 6 -10 7 
years. 

- to be surrounded by circumstellar dust (IR excess) and 
CO. CO detection indeed provides a good indication for 
the gas to be extended. 

- to have been observed at various wavelengths to allow 
comprehensive investigation. 

2.2. Observing strategy 

The high level of contrast between the investigated cir- 
cumstellar environments and their central star implies the 
use of coronagraphic techniques. The subtraction of resid- 
ual Point Spread Fonction (PSF) wings is required to take 
full benefit of the coronagraph's high contrast capabilities. 
For that, we observed in the same optical configuration a 
comparison star within a few degrees of the science target, 
of similar or slightly higher brightness, with similar spec- 
tral type and thought to be free of circumstellar matter. 

Comparison stars need to be observed close in time to 
each corresponding science target observation because of 
the slight possible variation of the PSF with time. This 
method has proved to be efficient for ground-based coro- 
nagraphic data with adaptive optics, especially because it 
is still efficient if the PSF is not perfectly symmetrical. It 
also allows the detection of circumstellar disks far from 
edge-on orientation. 

2.3. NICMOS2 data 

We obtained coronagraphic images of our four targets 
with the HST/NICMOS2 camera between February and 
November, 1998. All the observations were performed in 
Filter F160W (A c = 1.6pm, AA = 0.4pm). For the three 
targets presented in this paper, the observing log is sum- 
marized in Table [y. The individual integration times were 
determined so as to stay below the saturation limit of the 
detector on the brightest pixels (near the central hole edge 
area). Filter F160W was selected for optimum detectabil- 
ity: the PSF is well sampled and the background (zodiacal 
light and thermal background) remains low. 



2. HST/NICMOS2 observations and data analysis 

2.1. Target selection 

Our aim was to image circumstellar environments of tran- 
sition (old PMS) objects to establish the missing link be- 
tween disks around young embedded objects and (3 Pic-like 
stars in order to compare the disks with those (optically 
thick) of young stars on the one hand, and that (optically 
thin) of /3Pictoris on the other hand. Our four targets : 
HD 100546, SAO 206462, MWC 480 and HD 141569 were 
selected according to the following criteria : 
- to be PMS stars close to the MS, which can be character- 



2.4- Reduction procedure, photometry and astrometry 

For each observed star, the calibrated files (in counts/s) 
provided by the STScI are coadded to form a single image 
of the star during a same orbit. The reduction procedure 
then consists in subtracting the comparison star carefully 
scaled to the star of interest. The determination of the 
scaling factor is critical since a small change can signifi- 
cantly modify the photometry or in the worst case induce 
artifacts. It is assessed by azimuthally averaging the profile 
resulting from the division betwee n the star of interest b y 
the reference star images (see also Augereau et al. 1999b ). 
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Ta ble 1. Observing log. N ote th at SAO 206462 and S AO 206463 form a visua l binary system (HD 135344). References : 
<») |Waelkens et al. (1990cj jb|), m |Coulson ct al. (1995| ), < c > Simbad (CDS), < d > [Myers et al. (1998| ). 





Star Name 


Spectral Type 


V mag 


Date 


Integration 


Time 


Same Orbit ? 


Object 
Reference 


HD 100546 
HD 106797 


B9Vne 
A0V 


6.68 
6.07 


Nov 09, 1998 
Nov 09, 1998 


16x47.958 s = 
16x47.958 s = 


12m47s 
12m47s 


Yes 


Object 
Reference 


SAO 206462 
SAO 206463 


F4Ve<"» ; F8V< b) 
F2 (c) . A2<< i > 


8.65 
7.9 


Aug 22, 1998 
Aug 22, 1998 


7x111.931 s = 
7x191.961 s = 


13 m 03 s 
22 m 24 s 


Yes 


Object 
Reference 


MWC480 
HD 29646 


A2 

A2V 


7.72 
5.73 


Feb 24, 1998 
Feb 24, 1998 


11x39.953s = 
11x31. 959s = 


= 7m20s 
= 5 m 52 s 


Yes 



Table 2. Some kno wn astrophysical paramet ers o f the target sources pres ente d in this paper. Refer ence s : (o) |van dcr 
Anckej et al. (1997|), <"' |Sylvester et al. (1996|), '"' |Zuckerman et al. (1995|) , < d > |Coulson ct al. (1995| ), < e > punkin et af 



(1997|)j w |van den Ancker et al. (1998[ ), (a » |Mannings ct al. (1997b|) 



Star Name Distance [pc] Luminosity [L ] 



Mass [M e 



e) 



T cS [K] 



Age [Myr] 



Other Name 



HD 100546 

SAO 206462 

MWC480 



103±e w 
84(f) . 10 o<= 

131±^ (/) 



32.4i|; 5 5 (») 



32.4±* 3 4 3 <« 



2.4±0.1<° 



2.2±0.3 (/) 



O O + 0.1 (g) 

z " 3 -0.3 



10470<°) 
6250< d > ; 6660< e > 

8710<« 



> 10<° 



9 tc + 1 - 5 (/) 
z -°-0.9 



6 < 3 ) 



SAO 251457 
HD 135344 
HD 31648 



We use the 2.077x 10 -6 Jy . s / counts factor to convert 
NICMOS count rates to absolute fluxes and a zero point 
flux density of 1040.7 Jy to convert to magnitudes (NIC- 
MOS data handbook, version 4.0, dated december 1999). 
Point source photometry is obtained by integrating the 
total flux within a 0.5" radius circular aperture, then ap- 



plying a 1.15 correcting factor to compensate for the flux 
which fell out of this aperture (NICMOS photometry up- 
date web page). 



ISO observations also evidenced similarities between the 
10 /zm emiss ion of the circumstell a r dust and that of c omet 
Hale-Bopp (|Crovisier et al. 1997| ; |Malfait et al. 1998| ). 

As commonly observed for other HAeBe stars, 
HD 100546 exhibits photometric, polarization and spec- 
troscopic variab i lity (|Grady et al. 1996 , 1997 ; van den An- 



cker et al. 1998|; |Yudin fc Evans 1998| ; |Clarke et al. 1999 ; 
Vicira et al. 1999| ). The interpretation of s uch events in 



Th<i_ 
the rai ge 0.075"-0.076" in average 



particular in terms of star-grazing comets (Grady et al 



pixel scale slowly varies with time but stays within 1997|; |Vicira ct al. 1999|) remains ho wever uncertain ( |La- 



grange et al. 2000^ |Bcust et al. 2000J ) 



3. HD 100546 



3.1. A Werbici Be star close to the ZAMS 



The HD 100546 star has been intensively studied since it 



was iden tified as a member of the HAeBe group (Hu ct 
al. 1989|) . Indeed, HD 100546 is a B9Vne star showing a 



Another issue concerns the very distribution of the ma- 
terial surrounding HD 100546. The disk detection in scat- 
tered light b etween 40 (0.4") and 200 AU (2") by |Pantm 
et al. (2000 ) reveals a ~ 50° inclined disk almost in the 
SE-NW direction. At A = 1.3 mm, an extension as large 
as 15-20" in radius and close to detection limit is also 



strong infrared excess peaked at about 25 fim (IRAS) due 
to circumst ellar material and then fu lfills the criterions 
proposed by Waters fc Waelkens (1998J ) to identify HAeBe 
stars. Some astrophysical parameters are summarized in 
Table || 

The position of the star in the HR diagram indicates 
that HD 100546 is close to the ZAMS leading to an es- 



reported by Henning et al. (1998) in the same direction. 
From the theoretical point of view, it is unclear whether 
the circumstellar material arranges within a single disk. 
For instance, the presence of an additional envelope is 
proposed to reproduce the full spectral energy distribution 



(Henning et al. 1994) and to explain specific spectroscopic 
events (Vicira et al. 1999). 



timated age larger than 10 Myr (van den Ancker et al 
1997). This star, associated with the dark cloud DC 296.2- 
7.9 (Hu ct al. 1989|; |Vicira et al. 1999|), is s urrounded by 



a large amount of dust: between a few tens (Bouwman ct 



al. 200C) and a few hundred Earth masses (Henning et 



al. 1998). The dus t is rich in C and O and in particular 



PAH and silicates ( Malfait et al. 1998 ). The latter are very 
valuable since they constrain the dust optical properties. 



3.2. Results 

3.2.1. Resolved disk and vicinity 

The azimuthally averaged radial profile of the ratio be- 
tween HD 100546 and the reference star HD 106797 shows 
two basic distinct regimes: up to 3.5-3.8", the profile de- 
creases with distance from the star then reaches a plateau 
further out. The same behavior arises using the A2 and Al 
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[mog/Arcsec 2 ] 

Fig. 1. Reduced image of the HD 100546 disk observed in coronagraphic mode at 1.6 /im with the HST/NICMOS2 
camera. The coronagraphic mask is 0.7" in diameter (~ 9 pixels) whereas the numerical mask in the Figure is 0.9". 
Isophotes from 11 mag/Arcsec 2 to 17mag/Arcsec 2 are indicated. The dashed lines give the width and the position of 
the spider arms. The small panel in the upper right corner shows the result of isophotes ellipse fitting. Inhomogeneities 
due to the detector or to secondary spider diffraction spikes (which produce some kinds of fingers) are responsible for 
black spots. The most noticeable are labeled (a) to (i). In the same small panel, the plain lines represent the major 
and minor axis of the observed structure and the dashed ones represent the position of the spider arms. 



reference stars dedicated to MWC 480 and HD 141569 re- 
spectively. A resolved structure around HD 100546 is then 
detected below 3.5-3.8". A 3<r uncertainty of 2% is ob- 
tained on the scaling factor used to subtract the reference 
star to HD 100546 which directly impacts on the precision 
of photometric measurements. 

The final calibrated image is shown in Figure [P and 
reveals an extended and elliptical structure centered on 
the star. We interpret this structure in terms of circumstel- 
lar material located in an inclined disk. A visual inspection 
suggests the direction subtended by the spider arm at a 
position angle (PA, measured east of north) of ~ 148° 
as an axis of symmetry for the inclined disk, at least at 
large distances. A more robust method (ellipse-fitting to 
the isophotes) leads to a PA of 161°±5° but the result 
might be affected by different noise sources. Whatever the 
method, the global direction of the disk is consistent with 
Pantin et al. (20"00|) observations but the precise PA is 



assuming th at the disk is axisym metric, and agrees with 
the results of|Pantinet al. (2000| ). We measure a total flux 



found to be at leas t 20° larger than the measurements of 
Pantin et al. (2000| ). The axis ratio corresponds to a disk 
inclination of 51°±3° with respect to the line of sight, 



density of 73±7mJy (the uncertainty is dominated by the 
uncertainty on the scaling factor used to subtract the PSF 
star). 

Some companion stars are also detected in the close 
vicinity of HD 100546. The three brightest ones are 
shown in Figure [fl and measured astronomical parame- 
ters, namely projected separations, PAs and magnitudes, 
are summarized in Table H. 

3.2.2. Radial surface brightness profiles 

The radial brightness profiles of the disk are assessed by 
azimuthally averaging the surface brightness on elliptical 
contours. The radial profile measured towards the major 
axis of the disk is shown in the left panel of Figure 0. It 
reveals a smooth and continuous decrease of the surface 
brightness with the distance from the star from the very 
close edge of the coronagraphic mask up to 2.5". Beyond 
2.5", the profile is steeper and reaches the detection limit 
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Fig. 2. The HD 100546 disk radial surface brightness profiles azimuthally averaged on elliptical contours assuming 
PA = 161° for the major axis and an eccentricity of e = 0.63. The left panel shows the major axis profile and the right 
ones the minor axis profiles. The dark grey areas in the small boxes in t he leg end mimi c the angular sectors centered 
on the star used in each case. These results are discussed in paragraphs 3.2.2 and 3.2.5 . 



Table 3. Projected distances, position angles (PAs) and 
magnitudes in filter F160W (A c ~ 1.6 /im) of the 3 bright- 
est visual companion stars detected in the vicinity of 
HD 100546. The labels refer to Figure [|. 3a uncertainties 
on the positions and PAs are mainly due to the uncer- 
tainty on the central star position below the mask but 
also to the uncertainty on the companion star centroid. 
The 3<7 photometric uncertainty is induced by the uncer- 
tainty on the scaling factor used to subtract the PSF star 
to the star of interest (subsection 2.4). 



Label 


Projected distance 


PA 


magnitude 




[Arcsec] 


[degrees] 


filter F160W 


B 


4.50 ±0.15 


197.7 ±1.1 


15.74 ± 0.04 


C 


10.80 ±0.15 


277.1 ±0.5 


15.79 ±0.11 


D 


5.09 ±0.15 


156.3 ±0.9 


16.09 ±0.08 



(grey area in Figure ||) at 3.5-3.8". Basic radial power- 
laws properly fit this measured surface brightness profile 
with indexes -2.92 ± 0.04 in the radial range [0.5",2.5"] 
and -5.5 ± 0.2 outside of 2.7". 

The surface brightness profiles along the minor axis 
(NE and SW sides of the disk) are shown in the right 
panels of Figure 0. To allow direct comparison with the 
results along the major axis, we plot the profiles versus 
deprojected distances. They correspond to the observed 
distances along the minor axis over sin(51°). As a confir- 
mation of the isophotc ellipse fitting procedure, the main 
shape of the minor axis profiles superimposes well on the 
major axis one. 



3.2.3. NE-SW brightness asymmetry 

However, a significant brightness enhancement arises in 
the NE direction mainly below 0.7-0.8" in projected dis- 
tances but not in the opposite direction. Also in the same 
radial range, the behavior of the SW profile better follows 
the major axis profile down to ~0.5". Below this distance, 
the SW side of the disk is slightly fainter but with a poor 
level of confidence (very close to the detection limit). This 
NE-SW brightness asymmetry is clearly evidenced in Fig- 
ure |. 

To emphasize the asymmetry, we subtract a synthetic 
disk which fits the disk isophotes to the observed disk 
of HD 100546. In the NE side, the subtraction shown in 
Figure [| reveals an excess of flux up to 1.1 mag/ Arcsec 2 at 
0.5"-0.6" collimated in a direction close to the minor axis 
of the disk (plain line). In the SW side, the systematic 
pattern at about 2.2", responsible for the shape of the 
radial profile shown in the bottom right panel of Figure pi 
is clearly evidenced. 

We first checked whether this asymmetry could result 
from a bad centering of the PSF star during the reduction 
procedure. To significantly decrease the inner brightness 
asymmetry, the reference star has to be shifted by more 
than 1.5 pixel in the NE direction. This leads to unrealis- 
tic strong asymmetries in all the rest of the image which 
characterize an obvious bad centering of the subtracted 
PSF star. The asymmetry arises close to the position of a 
diffraction spike which might seem suspicious. But due to 
a good alignment, the diffraction spikes were completely 
removed during PSF subtraction and did not require any 
further process which would have directly cause the effect. 
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Fig. 3. Upper panel : Residue resulting from the subtrac- 
tion of a synthetic disk with an axisymmetrical surface 
brightness strictly identical to the one shown in the left 
panel of Figure g and inclined at 51° from edge-on to- 
wards PA=161° to the reduced image shown in Figure |l|. 
For each pattern, the lowest level is 0.4 mag/ Arcsec 2 . The 
white contours are spaced at 0.4 and 0.6 mag/ Arcsec 2 , the 
black ones at 0.8, 1 and 1.1 mag/ Arcsec 2 . As for Figure 
0, the plain lines represent the major and minor axis and 
the dashed ones the position of spider arms. Lower panel : 
radial profile along the NE minor axis averaged for each 
distance from the star over 9 pixels in the perpendicular 
direction. The subtraction highlights the NE-SW bright- 
ness asymmetry. It allows to quantify this asymmetry and 
to determine the PA of the excess. 



Could it be due to instrumental effects then? 
HD 100546 and the reference star have been observed dur- 
ing the same orbit. Therefore, orbit-to-orbit variations, 
such as the migration (~ 0.25 pixel) of the coronagraphic 



hole on the detector or defocus for instance (Schultz et al 
1999n7]cannot account for the observed effect. Another is- 
sue concerns the mis-centering of the star behind the mask 
during target acquisition. An accuracy of a few tenths of 
pixel is reached when centering PSF on the occulting mask 
but a mis-centering of the star by only a third of pixel can 
produce brightness asymmetries mainly within 0.4"-0.5" 
in radius (Schultz et al. 1999). Observations at different 



roll angles would help to check whether the asymmetry is 
real or an artifact. 



3.3. Discussion and interpretation 

3.3.1. Disk radiation to stellar luminosity ratio at 1.6 /im 

The measured flux density of the disk at 1.6 /im corre- 
sponds to 10.4±0.1mag according to the zero point flux 
provided by the S TScI. On the ot her hand, the V-H color 
index observed by Hu et al. (1989 ) and mor e recently Ku- 
rucz spectrum fitting by Malfait et al. (1998 ) indicate that 
the disk emission represents 55% to 60% of the total flux 
in the H-band. This excess may be due to the thermal 
emission of hot grains at very short distances from the 
central star, presumably close to the grain sublimation 
limit (a few fraction of AU, see also Figure ||). Obviously 
this excess falls below the coronagraphic mask used during 
present observations (the mask radius is 0.35" or ~36AU 
in situ according to Hipparcos measurements) . Assuming 
14=5.88 for HD 100546 ( [Malfait et al. 1998J ), our flux den- 
sity measurement provides a lower limit of 4 x 10~ 2 for 
the scattered to photospheric flux ratio at 1.6 /im with an 
uncertainty of 15%. This is about 2.5 less than the ratio 
inferred by Pantin et al. (2000|) in J band. 



3.3.2. Normal optical thickness of the disk 

As a first step in analyzing the data, we assume that 
the disk is optically thin whatever the direction. Since 
the disk is inclined with respect to the line-of-sight, it is 
then straightforward to derive the normal optical thick- 
ness rj_(r) at 1.6 /im. This requires some reasonable as- 
sumptions on the optical dust properties. Both method 
and assumptions are summarized in Appendix A. The re- 
sult is shown in the upper panel of Figure m In this sim- 
ple approach, the disk appears indeed optically thin ver- 
tically at 1.6 /xm everywhere in the resolved radial range 
but t± reaches a quite high maximum of ~ 0.06 close to 
the mask's edge. 
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Fig. 4. 1.6 (im normal optical depth derived from the 
observed radial surface brightness distribution of the 
HD 100546 disk in the optically thin approach. Between 
50 and 100 AU, it is a factor of between 40 and 50 less 



than the normal optical thickness inferred by Pantin et 
|al. (2000P at 1.25 /im (J band). 
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Following Krist ct al. (2000), we then assess the mid- 
plane optical depth by assuming that the vertical shape of 
the disk is a Gaussian with a similar scale height of 20 AU 
at 140 AU from the star and a (3 disk flaring index of 1. We 
find a 0.15 midplane optical depth between 50 AU and the 
outer edge of the resolved disk. It falls below 0.1 between 
~ 80 AU and the outer edge. This would indicate that the 
optically thin regime may not apply in the inner regions 
of the disk and that part of the stellar light is obscured at 
large distances in the midplane leading to an underesti- 
mation of t±_ . This is in lin e with the large millimeter flux 
( Hcnning et al. 1994 , 199S), that suggests that the disk is 
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at least marginally optically thick in part. 

3.3.3. Normal surface density and minimum mass of the 
dust disk 

The optically thin approach is usefull to determine a lower 
limit on the normal surface density a±(r) in solids (see 
Appendix B). Since we assume a grain size distribution 
proportional to a -35 , <Tj_(r) depends both on the mini- 
mum (a m i n ) and maximum (a max ) grain sizes in the disk. 
Dust grains are assumed to be icy porous aggregates made 
of a silicate core and coated by organic refractories. Two 
dust chemical compositions are considered: amorphous 
and crystalline grains with porosities of 0.5 and 0.95 re- 



100 200 300 400 

Distonce from the star [AU] 

Fig. 5. Normal surface densities of the disk which repro- 
duce, in the optically thin approach, the scattered light 
observations for different minimum grain sizes and the two 
types of grains considered. Note that surface densities as 
well as masses should be scaled by a factor y/a max [cm] if 
the maximum grain size was not 1 cm as assumed here. 



on the optical thickness in section 3.3.2, the inner part of 



spectively (see appendix A and Augereau et al. 1999a, for 
more details on the model used) . <r± (r) is plotted in Figure 
for different minimum grain sizes. Between 50 AU and 



100 At for instance, cr± (r) is as large as a few 10 b g.cm 2 



and about one order of magnitude larger if grains are 
at the blow-out size limit (~ 13 /im, see Figure A.l). In 



the disk can rather be considered of first generation. 

Outside of 270 AU, a±(r) falls as r - 3 - 5±a2 . A similar 
discontinuity of the dust distribution has already been re- 
ported in the case of the j3 Pictoris disk. If the dynamics 
of the HD 100546 disk further than 270 AU is similar to 
that of the outer part of j3 Pictoris disk (namely, supplied 
in small grains by collisions among planctesimals and ra- 
diation pressure), we would expect the disk radial density 
distribution to follow a —4 radial power law (Lecavclier des 
Etangs ct al. 1996J) . Such a disk, assumed optically thin 



and seen perfectly edge-on, indeed produces a midplane 
surface brightness proportional to r~ 5 (Nakano 1990) as 



terms of dust mass between 50 AU and 390 AU, this cor- 
responds to a few Earth masses at most (the dust mass 
is about 3750 x ax(50AU)M ffi ) and is at least one order 
of magnitude less than dust masses inferred from infrared 
or millim eter excesses ( Bouwman et al. 2000] ; Hcnning ct 



typically observed for (3 Pictoris (Heap et al. 2000). The 



latter conclusion also requires quite isotropic scattering 
properties. With these assumptions, the HD 100546 disk 



would then flare as 



„0.5±0.2 



But the dynamical similari- 



al. 1998|) . Obviously, if meter-sized bodies are formed, the 
normal surface density as well as the mass in solids should 
be multiplied by at least a factor of 10. The dust mass may 
also be significantly larger if the optical thickness inferred 
is underestimated. 



3.3.4. Shape of the dust distribution in the optically thin 
approach 

cr±(r) follows a —0.92 ±0.04 radial power law up to about 
270 AU. Similar distributions arc found for disks around 
TTauri stars (e.g. Dutrey et al. 1996). The modeling of 



cloud core infall with no viscosity also predicts disk surface 
densities proportional to r" 1 ( |Lin fc Pringlc 199C| ). Sim- 
ilarly, the surface density of the solar nebula during the 
disk dissipation is expected to follow a r _1 radial power 
law under the assumption that the utilization of available 
s olids to built p lanets in the Solar System is not efficient 
( |Cameron 1995| ). In that frame, and given the conclusions 



ties between the HD 100546 and (3 Pictoris disks at large 
distances remains an open issue. 

3.3.5. On the NE-SW brightness asymmetry 

What can produce the NE-SW brightness asymmetry if 
this is not an artifact due to instrumental effects? Al- 
though dust distribution asymmetries can not be ruled 
out, the fact that the brightness asymmetry of the 
HD 100546 disk occurs in a direction close to its minor 
axis tends to indicate that this would be rather due to 
scattering properties. 

In the case of an optically and a geometrically thin in- 
clined disk, anisotropic scattering properties can produce 
similar asymmetries but they should occur whatever the 
distance from the star. In present case, the asymmetry 
mostly appears at short distances where the normal opti- 
cal thickness reaches values close to 0.1. 

The transition between the (marginally) optically 
thick regime in the inner disk and the optically thin regime 
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at large distances combined with anisotropic scattering sponsible for most of the long wavelength IR excesses. 



properties may explain the shape of the observed NE-SW 
brightness asymmetry. Further observations and models 



At 100 /mi, the SED is mainly featureless and the IR ex- 



cess mostly comes from an underlying continuum (Mai 



are nee ded to confirm that issue. 



fait ct al. 1998). We can then compare the observed flux 



3.3.6. Dust temperature and minimum grain size in the 
disk 

Hu et al. (1989) first claimed that the HD 100546 dust 
disk consists of a cold and a hot dust populations to ac- 
count for the overall IR excesses. The relatively high spec- 
tral resolution of ISO spectra between A = 2.4 /mi and 



180 /mi yielded Malfait et al. (1998) to propose tempera- 
ture ranges for the two dust populations. Figure @ shows 
these four isotherms as a function of the distance from the 
star and of the minimum grain size a m i n in the disk. The 
light-grey areas represent then the two dust populations. 
We only consider amorphous grains since amorphous sili- 
cates are supposed to be at least ten times more numerous 



than crystalline ones (Malfait et al. 1998) 
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Fig. 6. Four mean isotherms of amorphous grains versus 
the distance from the star and the minimum grain size 
in the disk. The dark (resp. grey) isotherms correspond 
to the inner and outer boundaries of the cold (resp. hot) 
dust population proposed by Malfait ct al. (1998| ) from 
ISO spectra. The dashed area indicates the radial range 
in which the disk is detected in scattered light. 

In te rms of spatial exten t, the cooler population de- 
rived by Malfait et al. (1998| ) is consistent with the dust 
resolved in scattered light if the grains are larger than 
~ 0.5 /mi (see Figure |f intersection of the light-grey and 
dashed areas). Indeed, grains at T = 210 ± 10 K (the tem- 
perature corresponding to the inner edge of the colder 
population) smaller than a half micron lie further than 
~ 50 AU. Such grains would induce a detectable cut-off 
in the dust distribution somewhere between 50 AU and 
100 AU which is not observed. 

3.3.7. Consistency of the scattered light image with IR 
excesses 

The dust res olved in scattered lig ht, colder than ~ 210 K 
according to Malfait et al. (1998J ), is expected to be re- 



at 100 /nn to the thermal emission implied by the nor- 
mal surface densities deduced from present scattered light 
images in the optically thin approach (Figure 0). 
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Fig. 7. Contribution to the A = 100 /on thermal emission 
of the dust population resolved in scattered light down 
to 50 AU as a function of the minimum grain size a m i n 
(plain line). The dotted-dashed line corresponds to the 
same population but extrapolated down to 10 AU. The 
uncertainties on this estimation are: 1/ the ~ 10% uncer- 
tainty on the scattered image photometry, 2/ the uncer- 
tainty on the anisotropic scattering properties (we have 
assumed /(90°) = 1/4tt for simplicity, see Appendix A) 3/ 
all the molecules detected with ISO which could partially 
change the mean grain albedo and 4/ the exact surface 
density if the optical depth is underestimated. 

Figure R shows this comparison versus a m ; n for two 
different disk inner edges: the first one at 50 AU and the 
second at 10 AU. In the latter case, o\_ (r) has been extrap- 
olated down to 10 AU with a -0.92 radial power law. The 
distance of 10 AU has been chosen because it corresponds 
to the position of 210 K grains larger than typically the 
blow-out size limit (~ 13 /mi, Figures A.l and pft. 

Whatever a m i n , the predicted flux never reaches the 
observed limit. The optically thin approach does not then 
bring any upper limit on the minimum grain size in the 
disk. Even grains larger than 10-20 /jm lying in the 10 AU- 
inner-edge disk radiate a flux 2.5 to 3 times smaller than 
the measured emission. 

As already noticed, the optically thin approach may 
underestimate the surface density, indicating that part of 
the dust mass does not contribute to scatter light images 
but may emit at 100 /mi. For instance, with the assumed 
dust distribution, increasing the dust mass by a factor of 
40 would still be compatible with the observations if the 
grains are submicronic. 



4. SAO 206462 
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HST/NICMOS2 camera. 



4-1- The HD 135344 multiple system 



The F4Ve star SAO 206462 flWaelkcns ct al. 1990a| ]b] 



forms a visual binary system named HD 135344 with 
SAO 206463 -20.4" away almost in the Northern di- 
rection. The latest star, that we used as the compari- 
son star during our HST/NICMOS2 observ ations, has a 
spectral type between F2 (CDS) and A2 ( Myers ct al 
199g()~|Among our targets, SAO 206462 is the less lumi- 
nous star but exhibit s the largest relative infrared ex- 
cess (Lir/L* = 0.67, Coulson et al. 1995 ). Dust emis- 
sion features at 3.29 /im, 7.8 /zm a nd 11.3 /zm have been 
evidenced by Coulson et al. (1995 ) and attributed to aro 
matic molecules of a fe w angstroms ( Sylvester et al. 1996]) 
CO is also detec ted by Zuckcrman ct al. (1995 ) and Coul- 



son et 



al. (1998] ) whose abundance suggests a low gas to 
dust ratio (8 x 10~ 3 , Coulson et al. 1998| ) if CO is represen- 



tative of the gas content (i.e. not depleted with respect to 
the other species) consistent with that of Vega-like stars. 



Coulson et al. (1995) and [Malfait et al. (1998J ) derived 
from the full SED modeling a two dust component descrip- 
tion of the disk. However, without imaging, the cooler dust 
that is present in the outer part of the disk (> typically 
10 AU) is not precisely constrained. Moreover, as shown 
by Kenyon & Hartmann (1997), a single thin (not flared) 
disk can not account for the observed infrared excess and 
yields pylvester et al. (1997 ) to propose either the presence 
of an envelope or accretion activity. 

4-2. Results 

4.2.1. Disk detection limit 

No obvious feature that could be attributed to circumstcl- 
lar dust is detected in the 1.6 /im images of SAO 206462 
(Figure M). The corresponding detection limit (Figure 
H) shows that at 1.6 /zm the disk is fainter than 13.5- 
14mag.arcsec~ 2 further than 0.4" (~ 35 AU) and fainter 
than 16mag.arcsec~ 2 at 93 AU (1.1"), i.e. cl ose to the in- 



ner disk edge derived by Sylvester et al. (1997 ) from SED 
modeling. In the optically thin approach, this can be in- 
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Fig. 9. Upper limits on the disk normal optical thickness 
at A = 1.6 /zm (upper panel) and on the dust surface den- 
sity as a function of the minimum grain size a m i„ (lower 
panel) for SAO 206462 assuming that the disk is optically 
thin. 



terpreted in terms of an upper limit on the normal optical 
thickness of about 10~ 2 (Figure 0) outside of 0.4" in ra- 
dius (35 AU). The dust material further than 35 AU (if 
any) then is either optically thin at 1.6 /zm or does not 
contribute at this wavelength due starlight occultation by 
the dust closer to the star. 



4.2.2. Implications on the dust distribution 

With the rough assumption that the dust surface den- 
sity is proportional to r -1 ' 2 (as typically observed for 
young stars and which implies the larger mass), the dust 
mass between 35 AU and 420 AU would be less than: 

upper limit 



2270 x ct" 1 



35AU 
upper limit 



It corresponds to — O.25M0 assum- 



4 g.cm 3 . This must be compared to 



M 

mg ar;7=£Au ~ 10 
t he few Earth masse s required to reprod uce the full SED 
( Sylvester et al. 1997 ; Coulson et al. 1998 ) assuming more- 
over a maximum grain size at least 10 times smaller than 
for present work. 

A large fraction of the dust mass (optically thin ap- 
proach) or at least a large amount of dust (optically thick 
approach) is then certainly confined within the first few 
tens of AU. This is consistent with the mean cold popula- 
tion posit ion of 23 AU derived by Coulson ct al. (1995J ) but 
not with [Sylvester ct al. (1997 ) (who assume very small 
grains in the disk, — 50 A) . Present r esults agree with the 
disk size limit of — 60 AU inferred by layawardhana et al 
(2000D from their unresolved 10 and 18 /im images. 
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4.2.3. A tight binary system close to SAO 206462 

A close binary system in the SW direction is detected 
in the vicinity of SAO 206462 (Figure [u]). The bright- 
est companion is 5.8"±0.15" away from SAO 206462 at 
PA=129.8°±0.8°. The two close companion stars are sep- 
arated by only 0.32"±0.04" almost aligned along the NS 
axis; we find a PA of 352°±7 for the third compan- 
ion with respect to the brightest one. In Filter F160W, 
the secondary companion has magnitude of 15.2±0.1 and 
the third companion is about 6.3 times fainter (mag = 
17.2±0.1). 
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Fig. 11. MWC 480 detection limit 
HST/NICMOS2 camera. 



at 1.6 /an with the 




Fig. 10. The close binary system detected in the vicinity 
of SAO 206462 at 1.6 /xm with the HST/NICMOS2 cam- 
era. 



5. MWC 480 



5.1. A 



sua- arcs 



ec dust disk 



Among our sample, MWC 480 is the only one for which cir- 
cumstellar gas has been imaged ( Mannings et al. 1997b| ). 
The resolved structure shows Doppler shifts characteris- 
tic of the disk rotation over angular scale of a few arcsec 
and inclined by ~ 60° from edge-on. Meanwhile, the dust 
continuum emission in millimeter also appears elongated 
in the same direction but is less extended (HWHM ~ 0.6 " 
or -80AU). 

Two dust populations seem required to account for 
the overall infrared excesses. For instance, a combination 
of a geometrically thin disk and an envelope is proposed 



by Miroshnichcnko et al. (1999 ) (see also Malfait et al 



199S). As in the case of /3Pictoris, the 10 /im spectrum of 



MWC 480 shows features attributed to silicates (in par 
ticular crystalline olivine) similar to that of long-period 
comets in our solar system (e.g. comet Levy, [Bitko et al 
1999J). 



5.2. Detection limit and implications on the dust 
distribution 



According to the detection by Mannings et al. (1997b), 
the dust disk is expected inside an angular region smaller 
than 1" in radius, a region where it is challenging to resolve 
faint structures because of PSF subtraction residues and 
diffraction spikes for example. Actually, very close to the 
edge of the mask, we do not find any convincing features 
which would clearly reveal the presence of the circumstel- 
lar disk. The azimuthally averaged noise on the image, 
interpreted in terms of detection limit, is shown in Figure 

0. 

We applied in the case of MWC 480 strictly the same 
technique as the one used to derive some interesting prop- 
erties on the SAO 206462 disk. Similarly, we conclude that 
either the disk further than 60 AU is optically thin at 
1.6 /xm or that the starlight is at least partially masked 
by a large amount of dust inside 60 AU (Figure [lj) • 

For instance, if we constrain the dust mass between 
60 AU and 660 AU to be s maller than 170M e in s olids 
and assuming a^ ex r^ 1 - 75 ( Mannings et al. 1997b| ), we 
find an upper limit on the surface density at 60 AU of 
6 x 10~ 2 g.cm~ 2 . This is consistent with the results in 
the optically thin app roach from our scattered light ob- 
servations (Figure |l2|) even if grains are as large as a few 
hundred microns and even if we rather assume that only 
10% of the dust mass lies further 60 AU. The present ob- 
servations confirm then that most of the the dust mass 
is confined inside the fi rst hundred of AU as shown by 
Mannings et al. (1997b]) . 



6. Conclusion 

Among our four HST/NICMOS2 Cycle 7 targets, two cir- 
cumstellar disks have been resolved in scattered light. In 
this paper, we report the positive detection of an extended 
inclined disk around HD 100546 a typical transient old- 
PMS star expected to be in a disk dissipating phase. The 
shape of the derived surface density, the normal optical 
depth and, if confirmed, the brightness asymmetry along 
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Fig. 12. Upper limits on the MWC 480 disk optical thick- 
ness at A = 1.6 /xm and on the dust surface density as a 
function of the minimum grain size a m i n . 



the minor axis of the inclined disk indeed indicate that 
the HD 100546 circumstellar environment certainly traces 
a stage between "first generation" massive disks around 
PMS stars and those of "second generation" around MS 
stars (Vega-like disks). 

Disk non-detections are reported for both SAO 206462 
and MWC 480. In each case, the detection limit in scat- 
tered light allows to infer constraints on the dust location 
consistent with observations at other wavelengths. 

Acknowledgements. We thank the referee, Dr. M.L. Sitko, for 
helpful suggestions. 

Appendix A: Disk normal optical thickness 

The disk normal optical thickness can be straightforwardly 
derived from the normal surface brightness if we assume 
that the disk is axisymmetrical and optically thin. 

More precisely, the disk normal surface brightness in 
scattered light is: 

SB ± (r) ~ $ A r- 2 /(90°) x (a sca ) Q n ± (r) 

where: $a is the total received flux (star+disk) at the 
considered wavelength A, r is the distance from the star, 
/ is the scattering phase function, (er sca ) a is the scatter- 
ing cross-section averaged over the grain size distribution 
and nj_(r) is the grain normal surface density distribu- 
tion. This equation assumes that grain properties do not 
strongly depend on r. 

On the other hand, the normal optical thickness is 
given by: 
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Fig. A.l. Correcting factor e at A = 1.6 zmi as a func- 
tion of minimum grain size a m ; n (see text). The ratios of 
radiation pressure to gravitational forces /3 pr have been 
computed for grains with size a m i n in the vicinity of 
HD 100546, MWC 480 and SAO 206462 (assu med to be a 
F4V star with L»~2.5L Q and M« ~ 1.2 M , [Allen 1973|). 



with: 



(Cext/a 



where (u ex t) a is the averaged extinction cross-section. Fi- 
nally, we have: 



Ti(r) 



- 1 $3r 1 /(90 o )- 1 r 2 £B_L(r) 



The u> ratio depends on the wavelength, on the grain 
size distribution and on grain properties (shape, chem- 
ical composition ... ). In the following, we assume A = 
1.6 /im and a collisional grain size distribution propor- 
tional to a~ 35 between a min and a max . Grains are as- 
sumed to be spherical (a is the radius) and a max is fixed 
to 1 cm. We adopt the two types of grains as proposed 
in Augereau et al. (1999a), namely: amorphous "ISM- 
like grains" with porosity P of about 50% and crystalline 
"comet-like grains" with large porosity (P ~ 95%). Grains 
are assumed to be made of silicates, organic refracto- 
ries and a small amount of water ice (10%). We use the 
Maxwell-Garnett effective medium theory to compute the 
complex index of refraction of the aggregate and the Mie 
theory to derive the optical properties. 

For a m ; n ^> A, the lo ratio reaches a constant value ujq 
close to 0.5 ( Bohren fc Huffman 1983 ). Then if we note: 



T ± ,oW-2$- 1 /(90°)- 1 r 2 5 J B ± (r) 
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and 



: LOq X LU 



we have: 

T ± (r) ~ erj.,o(r) 

where e represents a correcting factor which must be ap- 
plied for the smallest a m i n , where the geometric optics 
approximation does not apply, e versus a mm is shown 



in Figure A.l for the two types of grains. Note that for 
a m in larger than a few 27 r(i~P) 1 / 3 ' w ^ s ver y c l° se to the 
albedo of the smallest grains. For simplicity, we assume 
that grains scatter isotropically, then: /(90°) ~ 1/An. 
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Appendix B: Disk normal surface density 

The normal surface density of the disk is obtained by: 

<r±(r) = (-irpa 3 ) a x n±(r) ~ C°"J-,o(r) 



with: <x±fl(r) = -p g T±. {r)Va min a n 
and: ( = ^ 



assuming luq ~ 0.5. The grain density p g depends on the 
grain properties. It is about 1.2g.cm~ 3 for amorphous 
grains and ten times smaller for crystalline grains due to 
the large porosity. 



The £ ratio plotted in Figure B.l is the correcting fac- 
tor (as e for the normal optical thickness). It must be ap- 
plied for the smallest a m i n , i.e. where the more numerous 
grains (those with sizes close to a m i n ) are not efficient scat- 
terers. But in the same time (for the same smallest a m i n ), 
£ is proportional to l/^/a m i n and then the normal surface 
density does not depend any more on the minimum grain 
size. 

The normal surface density aj_ (r) can therefore be in- 
ferred from the observed surface brightness if the disk is 
seen almost pole-on. If the disk is not resolved, it is also 
possible to derive an upper limit on uj_(r) from the de- 
tection limit. Obviously, the main uncertainty on a±(r) 
comes from the lack of strong constraints on the minimum 
and maximum grain sizes in the disk. 
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